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Abstract. - The Landau-de Gennes free energy is used to study theoretically the interaction 
of parallel cylindrical colloidal particles trapped at a nematic-isotropic interface. We find that 
the effective interaction potential is non-monotonic. The corresponding force-distance curves 
exhibit jumps and hysteresis upon approach/separation due to the creation/annihilation of 
topological defects. Minimization results suggest a simple empirical pair potential for the 
effective colloid-colloid interaction at the interface. We propose that the interface-mediated 
interaction can play an important role in self-organization and clustering of colloidal particles 
at such interfaces. 



Introduction. - Colloidal particles dissolved in a nematic solvent experience long-ranged 
interactions [1]. These effective interactions are generated by the distortions of the liquid 
crystal director around the particles and result in their clustering and self-organization [2, 
3]. Upon approaching and at the nematic-isotropic transition, the solvent forms nematic 
and isotropic domains, which leads to an additional interaction of the particles with the 
interface between the two phases [4,5]. The resultant morphology of the network formed by 
the particles in colloid-liquid crystal composites is very sensitive to this interaction [4,6]. In 
fact, colloidal particles can be captured by the nematic-isotropic interface and even dragged 
by a moving interface [4]. This enables one to manipulate tiny particles in suspension, which 
in a more general context is of technological importance for manufacturing of e-papers and 
electrophoretic displays [7] , separation of bacterial species and living cells [8] , trapping of DNA 
and polymer particles [9], as well as growth of photonic crystals [10]. The structures of self- 
assembled colloid layers at free surfaces of nematic films can be tuned by their thickness [11]. 

Recent experiments of dragging colloidal particles by nematic-isotropic interfaces [4] demon- 
strated that a single colloidal particle is attracted by the interface. This conclusion is in accor- 
dance with theoretical results obtained by minimizing the corresponding Landau-de Gennes 
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Fig. 1 - Aggregates of 16 /im polymer particles captured by a nematic-isotropic interface. The director 
is parallel to the nematic-isotropic interface and lies within the plane of the figure. In the experiment 
the interface moves in the upward direction and drags the colloids. Courtesy of A. Glushchenko, 
J. L. West, G. Liao, and Yu. Reznikov, Liquid Crystal Institute, Kent, Ohio, USA. 



Fig. 2 - Cross-section of two parallel cylindrical particles trapped at a nematic-isotropic interface at 
x = 0. The system is translationally invariant in z-direction. 



free energy [5] . The resulting force on the particle is non-monotonic and is roughly proportional 
to the particle radius. Experimental results [4] also suggested that the particles captured by 
the interface interact with each other differently than in the pure nematic or isotropic phase, 
e.g., the particle-particle interaction depends on the interfacial curvature; particles which in- 
teract only weakly in the bulk can segregate and form two-dimensional clusters once they are 
captured by the interface (see Fig. ^| . 

The appearance of additional interface induced features in the interaction between colloidal 
particles can be expected from analogous phenomena at interfaces of simple fluids [12] which, 
however, still pose unresolved challenges in understanding [13]. In liquid crystals the situation 
is even more complicated. The director deformations extend into the nematic phase, i. e., in 
addition to interfacial energies, the bulk elasticity contributes to the total free energy of 
the system. The effective surface tension (provided it can be introduced in the first place) 
depends on the orientation of the director, which in turn varies along the interface. The 
boundary conditions for the director at the particle surfaces give rise to topological defects 
accompanying the particles. In order to minimize the free energy associated with the director 
distortions as well as the free energy of the defect cores, the defects often merge with the 
isotropic phase [5] or with each other. 

The combination of the above-mentioned effects enriches the phenomena significantly but 
complicates their theoretical understanding enormously: expressions for the director field, 
the free energy, or force-distance profiles can hardly be obtained analytically. Numerical 
calculations are also not straightforward because the interfacial width and the size of the 
defect cores are normally much smaller that the size of the colloidal particle, so that there 
are several rather different length-scales involved. Therefore significant computational efforts 
involving finite-element methods with adaptive meshes are required to address even static 
problems [14]. 

In the following we study the interaction of two long parallel cylindrical colloids trapped 
at a (initially) flat nematic-isotropic interface. We find that their effective interaction is 
rather complex: the force-distance profiles can exhibit jumps and hysteresis upon approach 
or separation. The particles can either repel or attract each other depending on the material 
parameters and the geometry. Finally, we propose a simple effective pair potential which can 
be exploited to study clustering and self-organization of colloids at NI interfaces. 
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Landau-de Gennes free energy. - Two identical colloidal particles, each of which we take 
to be a long cylinder of radius R with the symmetry axis parallel to the z axis, are immersed 
into a nematic liquid crystal at a separation d, see Fig. [21 The nematic order parameter tensor 
at the boundaries is fixed in a such a way that without the colloids a flat nematic-isotropic 
interface is formed at x = 0. The order parameter at the top wall is fixed to zero; the absolute 
value of the order parameter at the bottom wall is fixed to the bulk order parameter of the 
nematic phase at two-phase coexistence. The director orientation at the bottom wall is also 
fixed and corresponds to the preferred director anchoring at the nematic-isotropic interface, 
i.e., parallel or perpendicular to the interface. The axes of the particles are positioned at 
x = x . In the course of the free energy minimization Xq adjusts itself such that that the 
a;— component of the force exerted on the particles vanishes. 

The system is characterized by the Landau-de Gennes free energy [15] 

/r /, | ^2 i 

dV a(T - T*)Q lj Q jl ~ bQijQjkQki + ri Q,jQj.)~ + — Q- ,.>Q + -^Qij^QikM ,(1) 

where summation over repeated indices is implied and the comma indicates the spatial deriva- 
tive. The positive constants a, b, c are assumed to be temperature independent, and T* is the 
supercooling temperature of the isotropic phase. The constants L\ and L 2 are related to the 
Frank-Oseen elastic constants K\\ = K 33 = 9Ql(Li + L 2 /2)/2 and K 22 = §Q\L\j2 and Qb 
is the bulk nematic order parameter. The sign of L 2 determines the preferred orientation 
of the director at the nematic-isotropic interface. L 2 > (L 2 < 0) favors planar (perpen- 
dicular) anchoring. We introduce the dimensionless temperature r = 2Aca(T — T*)/b 2 . The 
bulk nematic phase is stable for r < tni = 1 with a degree of orientational order given by 
Q b = 6(1 + y/1 - 8r/9)/8c; Q b (r > 1) = 0. 

We consider both rigid homeotropic and finite anchoring boundary conditions at the 
colloidal surfaces. In the latter case the surface free energy density is taken as [16] f s = 
±wTr(Q - S) where w is the anchoring coefficient and S a p = S (e a ep — 8 a p/2>) is the pre- 
ferred tensor order parameter at the colloid surface with e as the unit vector along the easy 
axis direction. 

For the constants entering the free energy density (Eq. JQ)) we use typical values for a ne- 
matic compound 5CB [17]: a = 0.044 x 10 6 J/m 3 K, b = 0.816 x 10 6 J/m 3 , c = 0.45 x 10 6 J/m 3 , 
L\ = 6xlO~ 12 J/m, T* = 307K. The choices for L 2 will be given below. The nematic-isotropic 
transition temperature for 5CB is Tni = 308.5 K. At coexistence the nematic coherence length 

1/2 

(i.e., the thickness of the nematic-isotropic interface) is £ = (24L lC /b 2 ) ' w 10 nm which sets 
the smallest length-scale of our description. 

The equilibrium distribution of the tensor order parameter Qij is obtained by minimizing 
the free energy functional Q numerically using finite elements with adaptive meshing. The 
area L x 2L is triangulated and the functions Qij are linearly interpolated within each triangle. 

Planar interfacial anchoring. For this case we chose the anisotropy of the elastic 
constants as L 2 /L\ = 2, favoring director alignment parallel to the NI interface and in the 
plane z = 0. We use rigid homeotropic boundary conditions at the colloidal surfaces. Upon 
minimizing the free energy we have identified two stable orientational configurations. At 
sufficiently large separations d both colloids are accompanied by a defect of strength 1/2 
(Fig-EIb)). If the particles move closer, the two defects merge and form a single, but rather 
extended, defect positioned between the colloids (Fig. Ela)). We have also observed three 
metastable configurations. Here either one of the defects, Fig. Etc), or even both defects, 
Fig. EUd) and (e), merge with the isotropic phase: the system is trying to reduce its free 
energy by annihilating the defects. However, the interface has to bend and almost wrap 
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Fig. 3 Fig. 4 

Fig. 3 - Order parameter and director maps for stable {a,b} and metastable {c,d,e} configurations. 
Particle radius R/£ — 5; elastic constants favor planar anchoring at the nematic-isotropic interface. 
Solid lines indicate the local director orientation. There is translational symmetry along the z axis 
and the director lies in the xy plane. 

Fig. 4 - Effective pair potential U (d) per unit length for planar director anchoring at the Nl interface 
and homeotropic anchoring at the colloid surface, for both stable and metastable configurations with 
R/£ — 5. The insets depict the corresponding configurations by iso-order-parameter lines. 

the particles in order to reach the defects. This of course increases the interfacial energy, 
which is roughly proportional to the area of the interface. For small separations, among the 
metastable configurations only the one shown in Fig- Et c ) is possible: due to strong anchoring 
at the particle surfaces the nematic phase forms a bridge between the particles, preventing 
the isotropic phase from reaching also the second defect. For larger separations, the isotropic 
phase can extend partially around the colloids and merge with the second defect. This can 
occur asymmetrically, as shown in Fig. IHId), or symmetrically, see Fig. Efe). 

The effective pair potential U(d) defined as the free energy per length of the colloidal 
particles is presented in Fig.0] In all cases the particles repel each other. The branch of U(d) 
which corresponds to the stable configurations, reveals practically no breaks in slope, even 
though at d/R « 0.4 the director field around the particles changes abruptly between the 
configurations (a) and (b). However, if one of the metastable configurations is present in the 
system, transitions between different branches can occur giving rise to jumps and hysteresis 
in the force-distance curves upon particle approach/separation. 

Homeotropic anchoring at the interface. - The choice L2/L1 = —1/2 for the anisotropy 
of the elastic constants favors director alignment perpendicular to the NI interface in the 
plane z = 0. In this case we have found only two, stable configurations. The first one, 
containing a defect, is shown in Fig.lSfa). It occurs at small separations, where the nematic 
bridge between the particles forces the director to be parallel to the interface. The other 
configuration has no defects, because the nematic bridge with director orientation parallel to 
the interface disappears and the defect merges with the isotropic phase. 

The corresponding effective potential U (d) is shown in Fig. El As in the previous case, 
the particles repel each other at small separations where there are large director distortions 
in the inter-particle region. However, contrary to the case with planar anchoring, there is a 
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Fig. 5 - Director and order parameter maps for homeotropic anchoring at the NI interface corre- 
sponding to elastic constants L2 = — l/2Li and particle radius R/£ — 5. Configuration (a) is stable 
at short distances, (b) at large distances. 

Fig. 6 - Effective pair potential U(d) per unit length for homeotropic anchoring at the NI interface, 
Z/2 = —1/21/1, homeotropic anchoring at the colloid surface, and R/£ = 5. The configuration with 
the defect results in strong attraction at short distances. There is repulsion at large distances. The 
insets characterize the corresponding configurations by iso-order-parameter lines (see Fig. [£). For 
d < d c configuration (b) is metastable. 



strong attraction between the particles for dj R < 1, where the configuration with the defect 
is stable. The branch without the defect results in weak repulsion, as for planar anchoring. 

Empirical expression for the effective potential. - Minimization of the free energy is a 
computationally expensive task, which cannot be used, for example, to calculate the forces on 
the particles each time their positions are advanced using molecular dynamics or Monte-Carlo 
algorithms for their kinetics. Use of effective potentials is a typical solution to this problem. 

In order to find such an appropriate potential we have calculated the free energy per unit 
length for particles of different radii, R/£ « 2.5 — 12.5, corresponding to R w 0.025 — 0.125 /Ltm. 
Since colloids are rather small, the finite director anchoring at the particle surfaces has to be 
taken into account. We have used a value appropriate for 5CB on a polymer surface [18], 
w£/Li = 0.822; for this energy the dimensionless anchoring parameter n = wQ\Rj K21 is of 
the order of 1. 

Minimization results for both planar and homeotropic anchoring of the director at the 
interface are shown in Fig. 0(a) and (b). The strong repulsion close to contact (d — 0) has 
become smoother for planar and disappeared for homeotropic anchoring at the interface, due 
to the finite director anchoring at the colloid surfaces. It turns out that in both cases the 
following shape of the potential fits the numerical data very well: 

r u a [i +1 {d/o 2 ], o<d<d c 

( )_ j [/(oo) + t/ (cxp [-a(d/£)~ /3 ] -l), d>d c , { > 

where a, /3, 7 are (positive) fitting parameters. The first part of the potential (0 < d < d c ) 
describes a simple harmonic-type attraction between the particles, which we observed for the 
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Fig. 7 - Effective pair potential as a function of the distance d between colloids for different colloid 
radii R: (a) planar and (b) homeotropic anchoring of the director at the interface and finite anchoring 
at the colloid surface. Symbols: minimization results; solid lines: fit to the empirical equation 
Inset illustrates the dependence of the fitting parameters on R. £ = 10 nm, K22 = 2.6 x 10~ 12 J/m. 
U a /K 2 2 ~ 6.8 - 26.5 corresponding to R/£ » 2.5 - 12.5. 



homeotropic anchoring; for planar interfacial anchoring d c = 0. In the case of homeotropic 
anchoring d c is a critical distance at which the defect annihilates (see Figs. Eland |BJ. 

We have no analytical support for the functional form of the effective potential |J2J. How- 
ever, qualitatively, the effective interaction of the colloidal particles in a bulk nematic phase 
is known to exhibit a power-law decay at large distances [1]. In our case, the leading term 
of the expansion for d/£ — * 00 is indeed ~ (d/£)~ /3 ■ On the other hand, if two colloids are 
sufficientlly close, so that the regions where the order parameter changes from their surface 
to the bulk value start to overlap, and where the defects (if present) interact with each other 
and the particles [14], the functional form of the potential changes. Normally, the distortions 
of the order parameter relax exponentially (similar to the variations of the order parameter 
through the nematic-isotropic interface) resulting in an exponentially varying strength of the 
interaction. This is also captured qualitatively by the potential given in Eq. J2J). 

The fitting curves are also shown in Fig. (a) and (b). For planar anchoring (3 rj 1.5; 
for homeotropic anchoring (3 « 1.0, i. e., in both cases the interaction potential is long- 
ranged. In case of homeotropic anchoring there is a critical size of colloidal particles (of the 
order of interfacial width) below which the discontinuous structural transition between the 
configurations with and without the defect disappears. Similar phenomena, named capillary 
bridging, can be observed in granular materials, liquid crystals, and binary mixtures [19]. In 
all cases, one of the phases (otherwise metastable in the bulk) condenses out and forms a 
bridge between the particles. This transition is discontinuous and ends in a critical point for 
small particle radii. In our case, the nematic phase located between the isotropic phase and 
the defect forms a bridge between the particles and results in their attraction. 

Conclusions. - In conclusion, within the Landau-de Gennes formalism we have studied 
the effective pair interaction between two cylindrical parallel colloids trapped at the nematic- 
isotropic interface. This pairwise interparticle force depends sensitively on the radii of colloids, 
the orientation of the director at the NI interface, and the anchoring strength, i.e., the nematic 
elasticity is an essential ingredient of the interaction. The formation or annihilation of defects 
influences the effective pair potential. It also changes the repulsive character of the interaction 
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at large separations to an attractive one at intermediate or small separations. Therefore, one 
should expect that many-body interactions are relevant to describe collective colloidal ordering 
at the interface. Another implication of our work is that the interface is able to rearrange 
the nucleation centers of, for example, phase-separating polymer network and help to design 
composite materials of extraordinary properties. 

* * * 

We thank to M. Deserno, A. Glushchenko, K. Kremer, G. Liao, M. Oettel, Yu. Reznikov, 
and J. L. West for useful discussions. 
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